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Abstract 

We investigate the influence of biquadratic exchange interactions on the low- 
lying excitations of a S = ^—ladder using perturbation theory, numerical di- 
agonalization of finite systems and exact results for ladders with matrix prod- 
uct ground states. We consider in particular the combination of biquadratic 
exchange interactions corresponding to ring exchange on the basic ladder pla- 
quette. We find that a moderate amount of ring exchange reduces the spin gap 
substantially and makes equal bilinear exchange on legs and rungs consistent 
with experimentally observed spectra. 
PACS numbers: 75.10Jm 
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I. INTRODUCTION 



Two-legged spin ladders with spins 1/2 have attracted considerable interest over the 
past years, both as ideal models for quasi ID materials and as a theoretical model for a spin 
liquid characterized by an excitation gap. From the geometric structure of the two-legged 
spin ladder as shown in fig. 1(a) it is clear that the main exchange interactions are expected 
along the legs and across the rungs; if the most general interaction for a plaquette formed by 
four spins S = \ on two neighboring rungs is considered, exchange along the two diagonals 
(corresponding to next nearest neighbor interactions) and biquadratic exchange interactions 
appear in addition. 

This general plaquette hamiltonian for S=l/2 is formulated as: 



For simplicity we have assumed equal exchange interactions J\ eg on the two legs and 
Jdiag on the two diagonals, a natural assumption for the symmetric structure of fig. 1 
which we will consider in the following. Recently theoretical studies have demonstrated 
that this generalized hamiltonian (i) allows to formulate models interpolating smoothly 
between the dimer and Haldane limits for the ground stated and (ii) has a parameter space 
sufficiently large to allow to study quantum phase transitions!!. The typical hamiltonian for 
both applications has some nonvanishing biquadratic terms. In the past, only little attention 
has been paid to biquadratic exchange terms, although finite strength of these terms was 
e.g. found in the spectra of small clusters of magnetic ions!. 

Recently, the possible importance of biquadratic exchange for some properties of low- 
dimensional spin systems has been pointed out: Honda et alE have argued for finite ring 
exchange (corresponding to a particular combination of biquadratic terms) on the basic CuC>2 




n 




) + Vll (S 




(1) 



2 



plaquette. The relevance of ring exchange for magnetization plateaus has been discussed for 
the spin ladder by Sakai and Hasegawai and for solid He 3 — films by Momoi et al.i. In this 
paper we study the effect of the biquadratic terms in the hamiltonian (1) on the low-lying 
excitations of this 'generalized spin ladder' by both analytical and numerical approaches. 
The motivation for this work is twofold: 

Firstly, we want to extend the knowledge obtained analytically in ref. ||] to a wider range 
of strengths of the biquadratic interactions. 

Secondly, we want to investigate the relevance of ring exchange of finite strength J T ing as 
introduced by Honda et al.i on the determination of coupling constants from ladder spectra. 
The standard analysis of experimental data on quasi ID ladder systems starts from eq. (1) 
with J diag and all biquadratic couplings set equal to and results in J ieg « 2 J rung from 
inelastic neutron scattering experiments on Sr 14 Cu240 4 ii and La 6 Ca 8 Cu 2 404ii, as well as 
from NMR experiments on these substancesi^S. This large value for the ratio Ji eg / J rung is 
not understood at present: It is not expected from the geometric structure of the ladder, from 
electronic structure calculations^ a value somewhat larger than unity is obtained and inter- 
ladder interactions^ which have recently been considered, do not resolve the discrepancy. 

In the following we show that even a moderate amount of ring exchange J r ing is relevant 
for the determination of the coupling constants in the two leg ladder from the energies of 
elementary excitations. We define the strength J r i ng on a plaquette of four spins i, j, k, I in 
terms of the permutation operator Pijki by the symmetric contribution 

Hring Jring ^ ] i^P'ijkl Plkji) 

plaquettes 

to the hamiltonian. Since we deal with spins S = 1/2, H r i ng is identical to a special choice 
of constants in eq. (1), namely 

Trring _ Turing _ rrring _ n j jring _ j jring _ jring _ 1 j 

V RR ~ v LL ~ V DD ~ z " J ring, J rung ~ J rmgi J\ eg — J diag ~ nn 9' 

A natural hamiltonian for the geometric structure of the spin ladder then is obtained by 
bilinear exchange terms and ring exchange terms, leading to the hamiltonian of eq. (1) with 
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°rung " rung ' ^ringi °leg leg ' ^^ringj " diag ^^ring 

VRR — VLL = = ^Jring- (2) 

In the following we will present results on the influence of J r i ng 7^ on the dispersion 
of low- lying excitations for antiferromagnetic ladders using perturbation theory (section II), 
numerical calculations covering the case of experimental interest (section III) and discussing 
the relevance for the dispersion of exactly known excitation branches (in the case of special 
choices for the strength of biquadratic exchange). 



II. PERTURBATION THEORY 

We have calculated the dispersion of the lowest triplet excitation in an expansion in 
the neighbourhood of the dimer point (where J rung = 1 is the only nonvanishing exchange 
constant)^. In the following we restrict the discussion to the hamiltonian of experimental 
interest 

l,n'S2,n + Jleg (Sl jn -Sl )n+ l + S2,n ' S2,n+1 ) } + Hring- (3) 
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The result for the spectrum of the lowest triplet excitation to third order in Ji eg , J ring is 
U)s=i(k) = (t + t\ cos k + t 2 cos 2k + £ 3 cos 3k) (4) 
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The leading coefficient for the dispersion is t\ oc J\ eg + J« ng in first order. This perturba- 
tive result shows: When J r ing > is present, an analysis of the spectra in terms of bilinear 
exchange only will lead to an effective value of Ji eg which is increased in comparison to the 
value found for J r i ng = 0. We will see in the next section from numerical diagonalization 



that this result of perturbation theory continues to be qualitatively true for larger values of 

Jieg- 



III. THE SYMMETRIC LADDER WITH BIQUADRATIC EXCHANGE 

Motivated by the results of recent inelastic neutron scattering experiments on the quasi 
ID ladder material LaeCagC^^Sr^i we have calculated numerically excitation spectra for 
spin ladders with the hamiltonian of eq. (2) using periodic boundary conditions for a finite 
number of spins. We have used the Lanczos method for ladders with a total of 24 spins i.e. 
12 rungs. Because of space inversion the number of wavevectors k with different energies is 
7, k = pn/6, p = 0, 1 . . . 6. In figs. 2a-d we present the results for the lowest negative parity 
excitation (parity with respect to interchange of the legs). This is the lowest excitation with 
the ladder gap at k = 7r; it is very likely a triplet excitation, since it is found with identical 
energies for both S* ot = 1 and S% ot = (S ta t > 1, however, is not formally excluded). We 
show the variation of the dispersion with J r ing between the limits —0.3 . . . + 0.3 for four 
values of Ji eg . It should be noted that the unit of energy is the bilinear rung exchange Jf ung 
(without the effective contribution J"no resulting from the ring exchange). J^ l ung should be 
fixed by the magnitude of the gap. 

For a simple quantitative presentation of the effect of ring exchange on the spectra we 
have considered the ratio 

R = ^. (5) 

This quantity may conveniently serve for comparison with experimental data since the 
most accurate spectra are obtained for k > | and it is here where uj(k) varies most. The 
curves in fig. 3 show the dependence of the ratio R on J ring for Ji eg / J rung = 0.5, 1.0, 1.5, 2.0. 
The analysis of the experiments gives R ~ 5, which may be realized for different combinations 
of Ji eg and Jring- Two examples are Ji eg ~ 1.8, J r i ng ~ and J\ eg ~ 1.0, J r ,i ng ~ 0.14. Thus 
we conclude from the data shown in fig. 2 the tendency that a finite J ring > tends to 
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simulate the effect of an increased bilinear exchange Ji eg on the legs. This is in agreement 
with perturbation theory as discussed in sect. II. We have checked that the perturbative 
result of eq. (3) is a good approximation to the numerical data for Ji eg = 0.5 for J r i ng < 0.1; 
for larger ring exchange a transition to a new phase occurs (see below) such that perturbation 
theory does no more apply. 

The main physics behind the importance of J r i ng for the determination of bilinear ex- 
change constants is the following: the influence of J r ing on the gap (at k = tt/2) is relatively 
much stronger than on the dispersion at intermediate wavevectors. Actually, for J\ eg < 1, 
the gap appears to vanish for a critical magnitude of Jring, indicating a phase transition. 
This phase transition may also be responsible for the irregular behaviour of the dispersion 
curves for Ji eg = 0.5, J r i ng = 0.3; this system may actually be in a different phase. It is the 
neighbourhood of this phase transition also for values J\ eg oc 0(1) which explains the strong 
influence of the biquadratic exchange terms. We expect this phase transition to be of the 
same type as discussed before for generalized spin laddersffll. leading to a spontaneously 
dimerized ground state. 

The energy of the spin gap is the basic input for the determination of the energy scale 
in ladder systems. Owing to the large effect of finite J r i ng on the spin gap the values of 
the exchange constants in Cu02— planes will have to be renormalized considerably. For 
J ring « 0.14 J rung , J leg « J Tung the spin gap is A « 0.28 J rung (compared to A « 0.51 J rung 
in the case J r i ng = 0). This implies an energy scale which is larger by a factor close to 2 
owing to the ring exchange terms. Thus the basic exchange energy J rung between two Cu 
ions changes from « 800K to ~ 1400/T, comparable to the magnitude of the main exchange 
constant in the 2D material La2Cu04. 

If the k— dependent dispersion is considered in more detail, the situation is of course more 
complex: In particular the presence or absence of a maximum between k = and k = n 
is affected by the value of J ring- At present, however, this is a point of little relevance for 
experimental results. In fig. 2b we have also included results for the lowest excitation energy 
with positive parity for J r ing = 0, ±0.3: these data illustrate our general observation that the 
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effect of biquadratic exchange on the positive parity excitation energies is less spectacular. 



IV. EXACT DISPERSION RELATIONS FOR SPIN LADDERS WITH SPECIAL 
VALUES OF THE EXCHANGE CONSTANTS 

It was observed in ref . [2] that the general plaquette hamiltonian for the ladder structure 
allows some combinations of parameters which result in either a rung dimer (singlet) ground 
state with exactly known triplet excitations or a Haldane-liquid like ground state with exactly 
known singlet excitations, both with negative parity. In both cases the dispersion is a pure 
cosine dispersion. The physical picture of these exact excited states is simple: An excited 
rung triplet (singlet) propagates to the neighboring rungs both on its right and left side 
without exciting states containing two or three quanta. The interesting and natural question, 
however, whether these exact excitations are of interest, in particular whether they are the 
lowest excitations in the corresponding ladder, can only be answered by numerical methods. 

The condition for an exact dimer ground state with exact triplet excitations for the 
hamiltonian of eq. (1) according to ref. |^ reads 




It has to be supplemented^ by the inequalities which guarantee the stability of the dimer 
ground state. For the hamiltonian of eq. (1,2) the conditions reads 

T bl - j ■ < - 

°leg ~ u ring SI g 

We illustrate this situation in fig. 4 where we plot the lowest negative parity excitation 
in the subspace S* ot = 1 for J\\ g = 0.15 and increasing values of J r i ng - We find that for 
J ring = Jii g = 0.15 the exact excitation is indeed the lowest excited state for Sf ot = as 
well as for ±1. The numerical spectra reproduce perfectly the exact dispersion law. Beyond 
the stability limit the exact excitation energy is between two highly excited states (since 
the rung dimer state then has an energy above the new ground state proportional to the 
number of rungs) and is therefore no more of interest. 
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It is interesting to note that an exact dimer ground state with an exact triplet excitation 
is also realized in 2D and 3D structures: The condition is that eq.(6) is fulfilled separately 
for the exchange interactions in each spatial direction. 

A second example where the exact ground state and exact excitation spectra can be 
explicitly given is the generalized Bose-Gayen model as introduced in section V of ref. |2|. 
This model is defined by the following choice of constants in the hamiltonian of eq. (1): 

J rung Vli Jleg 1; Jdiag 2/2 

Vrr = 0, V LL = ^(3 - 2y 2 ), V DD = ^(3y 2 - 2). (7) 

Varying the two parameters yi,y2 a number of phases, distinguished by different ground 
states are realized, among them the rung dimer and the valence bond (AKLT) phases. The 
phase diagram is reproduced in fig. 5. In ref. exact triplet (singlet) excitations in the rung 
dimer (AKLT) phase with negative parity have been given. We have investigated by exact 
diagonalization of ladders with 24 spins to what extent these exact excitations are the lowest 
ones for representative points in the phase diagram as indicated in fig. 5. Typical results of 
these calculations are shown in fig. 6. We find that close to the phase boundary of both the 
rung dimer (fig. 6a) and the AKLT (fig. 6b) phase the exact excitation is the lowest one for 
wavevectors in some neighbourhood of k = n (wavevector of the gap). This implies that they 
are exactly known critical modes as strongly suggested by physical intuition. Sufficiently 
far away in phase space from the transition lines, however, a crossover occurs and there 
may exist lower modes of both negative and positive parity (fig. 6c). On the line 2/2 = 1 
the exact excitations are dispersionless; this is reproduced for all wavevectors for the point 
(yi = 4,y 2 = 1), whereas for the point (yi = —4,y 2 = 0.95) the deviation from the exact 
solution is seen to occur for small wavevectors only. 
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V. CONCLUSIONS 



We have investigated biquadratic exchange interactions in their influence on the low- 
lying excitations of S — |— ladders with various coupling constants. From both analytical 
and numerical calculations we conclude that a quantum phase transition (which is likely to 
lead to a spontaneously dimerized ground state) occurs for small values of the ring exchange 
J ring > 0. At this phase transition the spin gap vanishes, therefore even small values of 
J ring imply large variations of the spin gap and have a strong influence on the determination 
of exchange parameters in the ladder in general. Quantitatively, the presence of a small 
amount of ring exchange is shown to be consistent with J rung ~ Ji eg as suggested by the 
geometrical structure of the ladder and with a value of this basic exchange constant between 
two neighbouring Cu ions which is close to the one found in the two-dimensional material 
La 2 Cu0 4 . 
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FIG. 1. (a) Structure of the two leg ladder, (b) possible ring exchange around the basic pla- 
quette. 
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J lBg =1 .0, both parities 




J le =1.5, parity=-1 





(d) 



G — GJ rlng =-o.3 
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J ri„ a = - 3 




FIG. 2. Ladder spectra (lowest negative parity triplet) for J r ing = —0.3 . . . + 0.3 and J\ eg = 
(a), 1.0 (b), 1.5 (c) and 2.0 (d). For J\ eg = 1.0 the lowest positive parity excitation (singlet) is 
given. 
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J, =0.15, parity=-1 




Jt/2 71 



k 

FIG. 4. Basic negative parity triplet for j\ l eg = 0.15 and J r ing = . . . 0.45. For J r ing = 0.15 the 
numerical results reproduce the exactly known dispersion curve. 
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FIG. 5. Phase diagram of the generalized Bose-Gayen model (see ref. 2) 
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y,=3.0 y 2 =-0.25 



y,=0.0 y 2 =0.27 




(c) 



y,=-4.0 y 2 =0.95 



— exact 
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- exact 
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FIG. 6. Numerical and exact excitations in the generalized Bose-Gayen model: The exactly- 
known triplet dispersions are numerically reproduced as lowest negative parity excitations for 
k > k m i n : (a): k m i n > 27r/3, (b): k m i n > 27r/3, (c): k m i n > tt/2. 



17 



